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PREFACE 

A feasibility  study  was  undertaken  at  the  U.  S.  Naval  Test  Pilot  School 
(USNTPS)  to  develop  a method  of  instrumenting  the  USNTPS  A-7C  airplanes  for 
flying  qualities  and  performance  testing  while  retaining  the  airplanes'  full  weapon 
systems  capability.  An  additional  purpose  of  this  study  was  to  minimize  the  cost 
A and  the  length  of  time  required  to  instrument  an  airplane  for  conventional  flying 

qualities  and  performance  testing.  This  technical  memorandum  presents  the  results 
of  this  feasibility  study.  The  work  was  conducted  as  an  integral  part  of  the  mission 
% of  USNTPS  under  NAVAIRSYSCOM  AIRTASK  A51 0-6 10C/053-F/6WSL7 1-0000  of 

2 December  1975. 
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INTRODUCTION 


BACKGROUND 


1.  A program  was  undertaken  at  the  U.  S.  Naval  Test  Pilot  School  (USNTPS)  to 
determine  the  feasibility  of  using  inertially  derived  data  to  generate  the  attitude, 
attitude  rate,  acceleration,  and  velocity  data  required  to  describe  aircraft  motion 
for  traditional  flying  qualities  and  performance  testing.  The  main  impetus  for 
pursuing  this  program  at  USNTPS  was  in  trying  to  retain  a full  weapons  system 
capability  in  an  A-7C  airplane  while  also  being  able  to  utilize  the  airplane  for  Test 
Pilot  School  flying  qualities  and  performance  testing.  Trying  to  incorporate  the 
currently  utilized  rate  gyros,  attitude  sensors,  accelerometers,  and  the  associated 
signal  conditioning,  etc.,  would  have  compromised  the  A-7C's  full  system 
capability,  and  the  cost,  in  terms  of  dollar  i and  aircraft  down  time,  would  have 
been  high. 


2.  For  many  years  aircraft  have  been  operating  with  highly  accurate  and  reliable 
inertial  navigation  systems.  When  one  looks  at  the  required  accuracies  of  the 
measured  accelerations  and  attitudes  from  an  inertial  system  that  navigates  with  a 
circular  error  probable  (CEP)  rate  of  1 to  2 nmi/hr  (1.9  to  3.7  km/hr),  one  finds 
acceleration  accuracies  of  mili  g's  and  attitude  accuracies  of  arc  minutes.  These 
accuracies  are  orders  of  magnitude  better  than  accuracies  with  conventional 
instrumentation. 

3.  This  program  was  not  the  f'rst  attempt  to  utilize  an  inertial  system  for 
aircraft  flight  testing.  The  General  Dynamics  YF-16  Flight  Test  Program  utilized 
a modified  Delco  Carousel  V inertial  platform  for  performance  testing  (reference 
1),  but  General  Dynamics  did  not  utilize  the  inertial  platform  for  stability  and 
control  testing.  This  platform  is  not  one  that  will  be  incorporated  in  the  production 
F-16  airplane.  General  Dynamics  was  highly  successful  in  this  venture  and  will 
again  utilize  this  inertial  system  for  the  performance  testing  of  the  preproduction 
F-16  airplanes.  Reference  1 discusses  the  use  of  an  inertial  system  for  the 
performance  testing  of  the  YF-16  and  should  be  required  reading  for  anyone 
intending  to  use  an  inertial  system  for  flight  testing. 

PURPOSE 


4.  The  purpose  of  this  feasibility  study  was  to  determine  if  an  inertial  navigation 
system  could  be  utilized  to  derive  the  conventional  flying  qualities  and 
performance  parameters  which  are  used  to  describe  an  aircraft's  motion. 

DESCRIPTION  OF  TEST  AIRCRAFT 

5.  A production  A-7C  airplane  was  used  for  the  feasibility  study,  and  it  was  an 
excellent  vehicle  with  which  to  pursue  an  inertially  derived  flying  qualities  and 
performance  parameters  (IFP)  program.  Even  more  important  was  the  availability 
of  an  A-7C  airplane  that  had  been  utilized  during  the  A-7C/E  weapons  system 
development.  This  airplane  had  the  capability  of  recording  or  telemetering  99 
parameters  from  its  on-  board  computer.  Included  in  these  parameters  were  those 
parameters  required  to  verify  the  feasibility  and  accuracy  of  an  IFP  system. 
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6.  The  A-7C  airplane  incorporates  a relatively  sophisticated  navigation  and 
weapons  delivery  system.  An  ASN-91  Tactical  Computer  (a  Four  Pi  IBM  solid  state 
general  purpose  digital  computer)  and  an  ASN-90  Inertial  Measurement  Unit  and 
Adaptive  Power  Supply  (built  by  Singer  Kearfott)  are  the  major  components  of  the 
airplane's  inertial  system. 

7.  One  ot  the  circuit  boards  (pages  is  the  IBM  term)  in  the  airplane's  tactical 
computer  was  modified  to  enable  the  recording  of  the  parameters  of  interest.  An 
intermediate  bandwidth  tape  recorder  to  record  a 50  kilo  bit  signal  was  also 
installed  in  the  airplane.  Since  this  initial  program  was  undertaken  only  to  prove 
the  concept  of  1FP,  no  additional  instrumentation  was  installed  in  the  airplane. 
Therefore,  cockpit  control  forces,  control  surface  positions,  and  engine 
performance  parameters  were  not  available  during  this  initial  feasibility  program. 

8.  The  A-7C  tactical  computer  has  a sample  rate  of  25  times  per  sec  for  most 
parameters  (all  inertial  parameters)  and  five  times  per  sec  for  others.  Table  I lists 
the  recorded  parameters  that  were  utilized  for  this  evaluation  and  their 
corresponding  sample  rates.  Of  the  99  parameters  from  the  tactical  computer  that 
are  recorded  on  one  track  of  the  tape  recorder,  only  18  were  utilized  in  this 
program. 
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Table  I 

Parameters  Recorded  in  Flight 


Parameter 

Sample  Rate 
(Times/Sec) 

Units 

Remarks 

Time 

25 

sec 

Identification 

5 

A one  (1)  appears  on  printout 
when  target  designate  feature 
is  activated;  serves  same  purpose 
as  "p"  light. 

Event  Number 

5 

' 

The  pilot  selects  the  desired 
event  number  by  selecting  the 
number  of  bombs  (0  to  99)  on 
Armament  Release  Panel. 

Inertial  Velocity  North  (V^) 

25 

fps 

Inertial  Velocity  East  (V_r) 
hi 

25 

fps 

Inertial  Velocity  Vertical  (Vyjj 

25 

fps 

Sine  of  Platform  Pitch 
Angle  (Sin  0 ) 

25 

- 

This  was  the  only  form  in  which 
pitch  was  available  from  the 
computer. 

Sine  of  Platform  Roll 
Angle  (Sin  -i ) 

25 

- 

This  was  the  only  form  in  which 
roll  was  available  from  the 
computer. 

Platform  Yaw  Angle  ( ; ) 

5 

deg 

! 

Vane  Angle  of  Attack  ( i) 

5 

' " 

deg 

i 

From  production  angle  of 
attack  vane. 

I 

Pressure  Altitude  (Hp) 

5 

ft 

Air  data  computer  (ADC) 
output. 

True  Airspeed  (V  ) 
A 

5 

fps 

1 

ADC  output. 

Smooth  Winds  East  (V™-,) 
Wls 

5 

fps 

Computed  by  tactical  computer. 

Smooth  Winds  North 

5 

fps 

Computed  by  tactical  computer. 

Magnetic  Heading 

5 

deg 

Mach  Number 

5 

- 

Air  Density 

5 

slug/ft3 

Weight  on  Gear 

5 

- 

3 


mm 
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9-  Th(>  ASN-90  is  a four  giinbaled  North  pointing  iner*ial  measurement  unit  which 
incorporates  integrating  accelerometers.  Previous  test  data  have  shown  the  system 
to  navigate  with  a mean  circular  error  probable  of  1.75  nmi/hr  (3  241  m/hr)  when 
operated  in  the  pure  inertial  mode  (reference  2).  The  inertial  system  incorporates 
a doppler  damping  feature  that  was  not  utilized  because  of  its  poorer  performance 
than  the  pure  inertial  mode. 

10.  A long-term  vertical  velocity  damping  feature  utilizing  the  air  data  computer 
unit  is  incorporated,  but  it  has  essentially  no  effect  on  short-term  inertial  vertical 
velocity  errors.  It  does,  however,  keep  the  d.  c.  or  long-term  value  of  vertical 
velocity  from  increasing  with  time. 

11.  The  IMU  is  rigidly  mounted  and  accurately  boresighted  to  the  body  axis  of  the 
airplane.  The  outer  roll  gimbal  of  the  IMU  is  shock  mounted  to  the  IMU  case  in 
order  to  attenuate  high  frequency  vibrations  and  to  provide  some  mechanical 
protection  to  the  unit.  The  IMU  accurately  measures  accelerations  up  to  a nominal 
input  fr  equency  of  20  cycles  per  sec. 

12.  Since  integrating  accelerometers  are  utilized  in  the  IMU,  a velocity  vice  an 
acceleration  output  is  obtained.  The  integrating  accelerometers  serve  a secondary 
function  of  attenuating  high  frequency  vibrations. 

SCOPE  OF  THE  TESTS 

13.  Two  flights  for  a total  of  3-3  flight  hours  were  flown.  On  the  first  flight  the 
inertial  navigation  system  (INS)  had  a failure  and  navigated  with  a 23  nmi/hr 
(43  000  m/hr)  error.  The  flight  was  therefore  repeated  after  the  INS  had  been 
repaired.  On  the  second  flight,  the  navigation  error  was  3 miles  (5  500  m)  after 
1.8  hours. 

RESULTS  OF  THE  STUDY 
ATTAINABLE  PARAMETERS 

14.  An  inertial  system  presents  attitudes  and  accelerations  (velocities  if 
integrating  accelerometers  are  utilized)  in  relation  to  its  three  orthogonal  axes 
which  are  maintained  stable  to  a desired  coordinate  system.  Typically,  this  is  an 
earth-referenced  "local  vertical"  system  having  one  axis  of  the  coordinate  system 
pointed  North,  one  East,  and  one  vertical. 

15.  In  flight  test  work,  an  earth-referenced  coordinate  system  is  utilized  when 
discussing  an  aircraft's  attitude.  However,  an  orthogonal  reference  system 
coinciding  with  the  aircraft's  body  axis  is  utilized  when  working  with  attitude  rates, 
accelerations,  velocities,  angle  of  attack,  and  sideslip.  Figure  1 shows  the 
earth -referenced  coordinete  system  and  the  aircraft  body-axis  coordinate  system. 
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16.  In  flight  test  work,  aircraft  motion  (a  vector)  is  described  in  terms  of  three 
components  in  the  orthogonal  triad  coinciding  with  the  aircraft's  body  axes.  Since 
the  inertial  system  presents  aircraft  motion  in  relation  to  an  earth-referenced 
orthogonal  triad,  it  is  necessary  to  transform  the  three  coordinates  of  the  aircraft 
■notion  vector  from  the  earth-referenced  frame  to  the  aircraft's  axis  referenced 
frame. 


17.  To  further  illustrate  the  need  for  this  transformation,  take  the  example  of  an 
airplane  executing  a turn  with  90  deg  of  right  bank  angle  and  its  nose  fixed  on  the 
horizon.  The  inertial  system  would  show  a platform  reference  pitch  rate  ( 0 ) of 
zero  and  a positive  platform  reference  right  yaw  rate  (ij').  However,  flight  test 
convention  would  describe  this  airplane  as  having  a positive  body  axis  pitch  rate  (q) 
and  a body  axis  yaw  rate  (r)  of  zero. 

18.  Euler  angles  were  utilized  to  transform  vectors  from  the  earth-referenced  axis 
system  to  the  aircraft's  body  axis  system.  Any  possible  rotational  sequence  could 
have  been  utilized,  but  since  aircraft  inertial  systems  have  their  gimbals  arranged 
from  inside  to  out  as  yaw,  pitch,  and  roll,  the  first  rotation  for  the  transformation 
was  about  the  platform  vertical  axis  (i.e.,  rotating  in  azimuth).  The  next  rotation 
was  in  platform  pitch  and  the  last  in  platform  roll.  It  is  important  to  always  utilize  > 
the  same  order  of  rotation  to  describe  the  attitude  of  the  aircraft.  With  a 
conventional  aircraft  INS,  this  will  always  be  yaw,  pitch,  then  roll. 

19-  Since  flight  test  work  references  aircraft  motion  to  the  air  mass, 
earth-referenced  velocities  must  be  converted  to  air  mass  reference  velocities. 
The  respective  North  and  East  component  of  the  wind  is  added  to  North  and  East 
inertial  velocities  prior  to  the  transformation  to  body  axis  velocities.  It  is  assumed 
that  the  wind  lay  only  in  the  horizontal  plane  and  that  the  vertical  component  of 
wind  is  equal  to  zero. 

AIR  MASS  REFERENCED  VELOCITIES  AND  ANGLES 


20.  Figure  2 presents  in  block  diagram  format  the  calculations  that  were  made  to 
the  INS  and  air  data  computer  (ADC)  outputs  in  order  to  obtain  the  air  mass 
referenced  velocities  along  the  aircraft's  three  orthogonal  body  axes.  The 
following  definitions  were  used: 


V„T,  Vy.  are  the  true  North,  East,  and  vertical  inertial  velocities 

obtained  from  the  INS. 


V V 
WN’  WE 

' Vt*  * ^ 


are  the  true  North  and  East  wind  velocities  as 
calculated  by  the  tactical  computer  by  comparing 
inertial  velocities  with  the  true  airspeed  from  the  ADC. 

are  the  earth-referenced  pitch,  roll,  and  yaw  angles  of 
the  aircraft  as  measured  by  the  INS  (i.e.,  platform 
pitch,  roll,  and  yaw  angles). 
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AIRCRAFT 

MOTION 


ADC 

TRUE  AIRSPEED 


Figure  2 

Transformation  Block  Diagrams 

21.  The  INS  inertial  velocities  (VNI,  V£I)  were  combined  with  the  tactical 

computer's  calculated  winds  (VWN>  VW£)  to  give  the  air  mass  referenced  velocities 

VN,  Vg,  and  Vy.  It  was  assumed  that  the  vertical  wind  component  was  zero. 

Therefore,  the  air  mass  referenced  vertical  velocity  (V^.)  was  assumed  equal  to  the 
inertial  vertical  velocity  (Vyj). 

22.  With  the  air  mass  referenced  velocities  computed  and  the  earth-referenced 
platform  attitudes  available  from  the  INS,  the  next  step  was  to  determine  the  air 
mass  velocities  relative  to  the  aircraft's  body  axis  coordinate  system.  The  matrix 
(reference  3)  used  for  this  transformation  is  given  below: 


V 

X 

cos  , cos  1) 

sin  y cos  0 

- sin 

< 

z 

V 

y 

= 

cos  sin  0 sin  y - sin  y cos  y 

sin , sin  • sin ; + cosy  cos  ! 

cos  sin  } 

V 

z 

cos  ysin  0 cos ^ + sin  14/ sin  (J 

sin  , sin 1 cos  T - cos  , sin  ; 

cos  1 cos ; 

> 

> 

The  solution  for  V , V , and  V is  then  as  follows: 
x y z 

V = V.,  cos  , cos  " + V,.  sin  C1  cos  • - V,.  sin  r 

x N E V 

V = VKI  (cos  , sin  sin  - sin,  cos:  } + V-,  (sin,  sin:'  sin,;  i cos  , cos  ; ) f V,,  cos  sin  : 

y N E V 

= Vj^j  (cos  sin  cos  ! + sin  , sin  : ) + Vg  (sin  . sin’.'  cos<  - cos  , sin  C ) + Vy  cos  cos  | 
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23.  The  next  series  of  computations,  illustrated  by  figure  3,  determined  the 
orientation  of  the  relative  wind  with  respect  to  the  aircraft  body  axis  coordinate 
system.  These  computations  were  utilized  to  determine  angle  of  attack,  a,  sideslip 
angle,  3 , and  flight  path  angle,  y . 


I.  d|  = tan 


' (vv») 

hd 

fe) 


I.  y = sin 


IV.  V 


AC 


= y/v  2 + V 2 + V 2 
V x y z 


Figure  3 

a , 3 » Y Determination  Block  Diagrams 


24.  During  dynamic  maneuvers,  rapid  accelerations/decelerations,  transonics, 

sideslips,  and  normal  accelerations  other  than  one,  static  source  position  error  and 

pitot  static  system  lags  can  cause  appreciable  errors  in  the  true  airspeed  derived  by 

the  ADC.  It  is  for  this  reason  that  the  computed  winds  were  "frozen"  in  the 

beginning  of  most  data  gathering  maneuvers,  and  a "corrected  true  airspeed"  (V . _) 

was  calculated  for  the  computation  of  3 and  y . The  mean  of  the  wind  readouts 

during  the  previous  5 sec  prior  to  freezing  the  wind  was  utilized  as  the  magnitude 

and  direction  of  the  frozen  wind.  Since  the  frozen  value  of  the  wind  was  used  for 

the  calculation  of  V and  V , the  computed  angle  of  attack  was  not  subjected  to 
x y 

the  pitot  static  measurement  errors  which  occur  during  the  various  maneuvers. 
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25.  The  calculation  of  winds  by  the  A-7's  on-board  computer  assumes  that  the 
aircraft  is  in  zero  sideslip  (i.eM  the  true  airspeed  vector  is  assumed  as  being 
attached  to  the  aircraft's  nose).  It  was  therefore  necessary  to  trim  the  airplane  for 
zero  sideslip  at  the  beginning  of  each  maneuver  prior  to  freezing  the  winds.  In  a 
symmetrically  loaded  airplane  with  good  side  force  characteristics  (d<J/dr),  zero 
sideslip  is  quite  accurately  obtained  by  trimming  the  ball  of  the  needle-ball 
instrument  to  the  center,  i.e.,  zeroing  any  lateral  accelerations  on  the  aircraft.  If, 
however,  the  airplane  has  weak  side  force  characteristics  or  has  an  asymmetric 
load  or  asymmetric  thrust,  an  independent  reference  to  give  zero  sideslip  is 
required.  A simple  yaw  string  could  be  utilized  and  would  be  expected  to  put  the 
airplane  within  ' 1 deg  of  zero  sideslip.  As  +he  subsequent  error  analysis  shows,  the 
error  in  initially  obtaining  zero  sideslip  acts  much  like  a tare  correction,  and  the 
calculated  sideslip  for  most  maneuvers  will  maintain  the  constant  initial  trim  error 
in  3 . For  the  case  where  heading  is  maintained  essentially  constant,  the  3 error 
will  be  equal  to  the  initial  zero  3 error  times  the  cosine  of  the  bank  angle,  and 
angle  of  attack  will  be  in  error  by  the  magnitude  of  the  initial  zero  3 error  times 
the  sine  of  the  bank  angle. 

26.  Of  course,  a conventional  sideslip  vane  on  a nose-mounted  boom  could  be 
utilized  to  give  zero  3.  However,  the  use  of  IFP  will  obviate  the  need  for  boom 
mounted  u and  3 vanes  if  another  type  of  zero  sideslip  reference  is  utilized  such  as 
a yaw  string. 

BODY  AXIS  ACCELERATIONS 

27.  The  acceleration  of  the  aircraft  with  respect  to  the  three  aircraft  body  axes 
(Ax>  Ay  Az)  was  determined  by  differentiating  the  inertial  velocities  (V^,  Vg,  Vy) 

and  then  transforming  the  resulting  accelerations  to  the  body  referenced  axis 
system.  A "flat  earth'’  assumption  was  made,  and  no  attempt  was  made  to 
compensate  for  the  inertial  system's  torquing  of  its  platform  to  keep  itself  aligned 
to  true  North  and  an  earth-referenced  vertical  and  to  compensate  for  the  coriolis 
effect.  Because  of  the  small  magnitude  of  the  resulting  errors  which  occur  at  the 
airspeeds  flown  by  a transonic  aircraft,  they  can  be  ignored  with  negligible  effect 
on  the  data.  A more  extensive  discussion  of  these  errors  is  presented  in  paragraphs 
55  through  77.  If  desired,  corrections  for  these  effects  can  easily  be  made. 

28.  Inertial  velocities  were  smoothed  over  five  data  bits  using  a third  order  least 
squares  curve  fitting  routine.  The  derivative  at  the  data  point  was  then  taken  in 
order  to  calculate  acceleration.  In  other  words,  the  two  velocity  data  bits  prior  to 
and  the  two  after  the  data  point  were  utilized  for  taking  the  derivative,  i.e.,  a time 
span  of  0.16  sec  when  data  were  computed  at  25  times  per  sec.  It  was  important  to 
minimize  the  time  span  over  which  the  data  were  smoothed  to  improve  data 
accuracy  during  dynamic  maneuvers. 
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FLIGHT  PATH  REFERENCED  ACCELERATIONS 

29*  Acceleration  along  the  aircraft's  flight  path  was  calculated  by  the  following 
equation: 

Acceleration  Along  the  Flight  Path  (AFP): 

AFP  = A cos  (X  + A sin  x + A sin  (•> 
x z y 

30.  In  order  to  include  the  earth's  gravity  component  and  have  parameters  that  are 
analogous  to  the  outputs  of  accelerometers  maintained  normal  to  and  along  the 
flight  path,  the  following  equations  were  used: 

Acceleration  normal  to  the  flight  path  (NZFP)  - Units  of  g's: 


A sin  v.  A cos  a 

NZFP  = — — - + cos  o cos  * 

g g 


Indicated  Acceleration  Along  the  Flight  Path  (AFPI)  - Units  of  g's: 


AFPI  = + sin  Y 

g 

31.  Having  less  flight  test  value  but  useful  for  comparing  IFP  data  with  that  from 
an  accelerometer  mounted  along  the  aircraft's  z axis  is  the  conventional  normal 
acceleration  (N  ): 

2 A 

N = - z + cos  0 cos  '{> 
z — 
g 

BODY  AXIS  ATTITUDE  RATES  AND  ACCELERATIONS 

32.  The  aircraft's  attitude  rates  and  attitude  accelerations  with  respect  to  the 
aircraft  body  axis  were  determined  by  the  following  equations: 
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33.  0,4),  and  ij,  were  determined  using  the  same  third  order  curve 
fitting/differentiation  routine  utilized  for  determining  body  axes  accelerations.  In 
the  same  manner  after  p,  q,  and  r were  calculated,  their  derivatives  were  taken  to 
yield  pitch,  roll,  and  yaw  accelerations  p,  q,  and  r. 

INERTIAL  MEASURE*  1ENT  UNIT  LOCATION 

34.  If  the  inertial  measurement  unit  (IMU)  and  C.  G.  are  not  collocated,  the  IMU 
will  experience  different  accelerations  than  the  C.  G.  of  the  aircraft.  Fortunately, 
the  IMU  in  the  A-7C  airplane  is  located  within  6 in.  (.15  m)  of  the  C.  G.,  and  C.  G. 
travel  with  fuel  usage  is  minimal.  Therefore,  no  compensation  was  made  for  the 
relative  location  of  the  IMU  with  respect  to  the  C.  G. 

35.  If  the  IMU  were  located  at  a position  other  than  the  C.  G.,  the  motion 
measured  by  the  IMU  can  be  transformed  to  motion  about  the  C.  G.  The  relative 
locations  along  the  x,  y,  and  z axes  between  the  IMU  and  C.  G.  must  be  known. 
Appendix  B presents  the  equations  which  can  be  utilized  to  transform  IMU  motion 
to  aircraft  C.  G.  motion,  or  in  fact,  the  motion  of  any  point  on  the  airplane. 

ACCURACY  OF  THE  DATA 

SOURCES  OF  ERRORS 

36.  The  objectives  of  this  section  are  to  identify  the  primary  IFP  error  sources  and 
to  show  how  each  of  these  sources  contributes  to  an  error  in  the  IFP  data.  The 
error  sources  in  the  IFP  data  can  be  categorized  into  four  groups: 

a.  INS  errors. 

b.  Air  mass  measurement  errors. 

c.  Flat  earth  assumption  errors. 

d.  Differentiation  errors. 

INS  Errors 


37.  The  technology  of  measuring  the  numerous  errors  of  interest  in  inertial 
systems  covers  a wide  range  of  specialty  fields,  hi  general,  there  are  two  opposing 
methods  of  measuring  these  errors.  One  is  to  make  a sufficient  number  of 
measurements  to  statistically  define  the  errors;  the  other  is  to  use  instrumentation 
with  accuracy  an  order  of  magnitude  better  than  the  errors  being  measured.  This 
last  requirement  challenges  the  accuracy  of  measurement  standards  (reference  4). 
Instead  of  actually  attempting  to  measure  the  INS  errors,  the  approach  taken  in 
this  section  is  to  show  how  a given  error  in  an  INS  (attitude  or  acceleration)  affects 
the  INS  output  (velocity  and  position). 
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38.  Two  examples  will  be  used  to  illustrate  the  relationship  between  the  INS 
attitude  and  acceleration  errors  and  the  resultant  navigation  errors.  It  should  be 
noted  that  the  A-7C's  INS  navigates  with  a mean  navigation  error  of  less  than  2 
nmi/hr  (3  700  m/hr). 

EXAMPLE  1 

39.  A single  axis  INS  is  modeled  as  an  accelerometer  on  a stabilized  level  platform 
(figure  4). 


VERTICAL 


NORTH 


40.  The  above  example  is  representative  of  the  North  axis  of  an  INS  with  the 
Schuler  feedback  loop  removed.  It  is  assumed  that  the  INS  initial  velocity  is  zero 
and  the  actual  North  acceleration  of  the  INS  is  zero  (AN  = 0).  However,  the 
accelerometer  is  assumed  to  have  an  acceleration  measurement  error  of  .001  g 
(AA^  = .001  g).  This  accelerometer  error  will  be  interpreted  erroneously  as  a 
true  acceleration,  and  the  INS  will  integrate  the  error  as  illustrated  below. 


aN 


1 — ; 

1/S 

L*J 

1/S 

V 

'N 


VN  - INS  VELOCITY 
PN  - INS  POSITION 


41.  Figure  5 presents  the  INS  velocity  error  (AV^)  and  position  error  (APjj)  as  the 
result  of  this  .001  g error  in  acceleration  measurement. 
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42.  In  the  above  example  a .001  g accelerometer  error  caused  velocity  and 
navigation  errors  of  116  fps  (35.4  m/s)  and  34  nmi  (63  km),  respectively,  after  1 
hour.  As  will  be  shown  subsequently,  the  short-term  velocity  error  is  of  primary 
importance.  The  maximum  velocity  error  chnge  in  a 2 minute  period  was  only 
3.86  fps  (1.18  m/s)  (or  3.3  percent  of  the  total  velocity  error  after  1 hour)  despite 
the  very  large  navigation  error.  If  ar  attitude  error  occurred  and  the  platform  was 
not  level  (i.e.,  6 / 0),  the  accelerometer  would  indicate  an  erroneous  acceleration. 
An  attitude  error  of  +.0573  deg  would  cause  the  accelerometer  to  indicate  an 
erroneous  .001  g.  The  INS  will  interpret  this  signal  as  true  acceleration  to  the 
North  and  consequently  produce  the  same  velocity  and  position  errors  as  depicted 
in  figure  5. 

EXAMPLE  2 


43.  A more  accurate  model  of  the  INS  which  includes  the  Schuler  loop  is 
illustrated  below  (figure  6). 


■VN 


VN  INS  VELOCITY 

PN  INS  POSITION 

R EARTH'S  RADIUS 

g GRAVITATIONAL  CONSTANT 
(32.2  FT/SEC2  or  9.81  M/SEC2) 


Figure  6 

(Adapted  from  Reference  3) 
Simplified  INS  Block  Diagram 


44.  Again  it  is  assumed  that  the  INS  initial  velocity  is  zero  and  the  actual 
acceleration  of  the  INS  is  zero.  A .0C1  g error  in  the  accelerometer  (acting  at 
point  A)  or  equivalently  a .057  deg  tilt  error  of  the  platform  (point  B)  would  cause 
identical  INS  velocity  and  position  errors.  Figure  7 presents  the  INS  velocity  and 
position  errors  resulting  from  the  .001  g acceleration  or  .057  deg  attitude  errors. 
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46.  The  above  examples  illustrate  that  the  navigation  errors  are  very  sensitive  to 
acceleration  and/or  attitude  errors.  Those  two  simplified  examples  of  an  INS 
should  show  that  an  INS  that  navigates  with  satisfactory  accuracy  (less  than  Z 
nmi/hr  or  3.7  km/hr)  measures  accelerations  to  accuracies  of  mili  g's  and  attitudes 
to  accuracies  of  arc  minutes. 

Air  Mass  Measurement  Errors 


47.  The  basic  vector  equation  representing  the  navigation  solution  is 


Where  - velocity  vector  of  the 
aircraft  relative  to  the 
air  mass 


Vy  - wind  velocity  vector 

(velocity  of  the  air  mass 
relative  to  the  inertial 
frame) 

Vj  - velocity  vector  of  the 
aircraft  relative  to  the 
inertial  frame. 


Solving  the  above  equation  for  V^.  gives 

= Vj  - Equation  A 

48.  The  above  equation  shows  that  the  wind  vector  can  be  obtained  by  taking  the 
difference  between  the  air  mass  reference  velocity  (V^)  and  the  inertial  velocity 

(Vj).  This  equation  is  mechanized  in  the  A-7C's  tactical  computer.  In  order  to 
investigate  the  effect^  of  errors  in  and  Vj  on  the  IFP  data,  the  "pleasured"  air 
mass  velocity  vector  (V^)  and  the  measured  inertial  velocity  vector  (Vj)  are  defined 


and 


V = V + A V 
I I U I 


Equation  B 


+ A V. 
A 


49.  The  above  equations  illustrate  that  the  measured  velocities  (Vj  and  V^)  contain 

error  vectors  (AVj,  AV^).  Therefore,  the  wind  vector  computed  from  the 

"measured"  inertial  velocity  and  air  mass  velocity  will  also  contain  errors.  That  is, 
the  computed  or  "measured"  wind  (V^)  is 

?w  = (Vj  + AVj)  - (V  - A VA)  Equation  C 
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50.  A comparison  of  the  measured  wind  with  the  ideal  wind  equation  (Equation  A) 
shows  that  the  wind  error  is 


AVW  = AVAVA 

51.  For  flight  test  work  one  is  usually  interested  in  determining  an  air  mass 
referenced  velocity  (V^).  The  normal  convention  is  to  express  the  magnitude  of 
the  vector  as  the  true  airspeed  (V^).  Hie  direction  of  the  vector  relative  to  the 
body  axis  system  is  specified  by  the  angle  of  attack  (a)  and  the  sideslip  (o).  The 
tollowing  development  shows  that  the  air  mass  velocity  vector,  computed  from  an 
inertial  velocity  vector  and  wind  velocity  vector,  is  not  dependent  on  the  constant 
portion  of  the  inertial  velocity  error. 


52.  The  wind  vector  was  computed  at  the  beginning  of  each  maneuver  with  the 
airplane  in  stabilized  flight.  The  computed  air  mass  velocity  is  thus 


A I 


W 


Substituting  Equations  B and  C into  the  above  gives 


V 


w 


Vj  + A V, 


I-[<VAVI> 


(vA  - ava) 


53.  Since  the  wind  is  frozen  at  the  beginning  of  many  maneuvers  (time  = t ),  the 
air  mass  referenced  velocity  is  subsequently  computed  from  the  measured  inertial 
velocity  and  the  "frozen"  wind.  Under  these  conditions,  the  "measured"  air  mass 
velocity  at  any  time  during  the  maneuver  becomes 


W(t  ) 
o 


.A- 


7,  ♦ AV, 


7i(.  ) * 'V  ) - (9A(.  ) 
o o o 


AVA(.  )> 
o 


54.  An  analysis  of  the  above  equation  yn  . the  following: 

a.  The  inertial  velocity  error  can  be  broken  into  two  parts:  a constant  error 
portion  and  a time  dependent  eiror  portion. 

1.  Computed  air  mass  refererced  velocity  is  unaffected  by  the  constant 
portion  of  the  inertial  velocity  errors  (i.e.,  when  AVj  = AVj^  ^). 


2.  Computed  air  mass  refe’enced  velocity  is  affected  by  changes  in 
inertial  velocity  error  during  the  maneuver.  Therefore,  A Vj 

A VKt J will  introduce  errors  into  the  computation  of  air  mass 

o 

referenced  velocity  (V^).  However,  this  term,  the  short  term  inertial 
velocity  error,  was  shown  previously  to  be  a very  small  error  source. 


b.  The  errors  in  the  air  mass  velocity  vector  can  be  broken  into  three  parts:  a 
wind  change  during  th<*  maneuver,  an  ADC  true  airspeed  measurement 
error,  and  the  airplane  not  being  in  zero  sideslip  when  the  winds  are  frozen. 
These  affect  computed  air  mass  referenced  velocity  as  follows: 


TM-TPS  76-1 


1.  Any  change  in  wind  magnitude  and  direction  during  a maneuver  when 
the  winds  are  frozen  will  cause  the  computed  air  mass  referenced 
velocities  (V  and  V ) to  be  in  error  according  to  the  corresponding 
magnitude  and  direction  of  the  wind  change. 

2.  If  the  ADC  incorrectly  computes  true  airspeed,  the  computed  wind  will 
be  in  error  by  the  magnitude  of  the  ADC  error  (Note  1).  The  direction 
of  this  wind  error  will  be  equal  to  the  reciprocal  of  the  heading  of  the 
airplane. 

3.  If  the  airplane  is  not  in  zero  sideslip  when  the  winds  are  frozen,  the 
wind  computation  will  be  in  error.  The  wind  error  will  have  a 
magnitude  equal  to  the  sine  of  the  sideslip  angle  multiplied  by  the  true 
airspeed  (V  sin  ;•' ) and  a direction  perpendicular  to  the  aircraft's  true 
heading. 

NOTE  1.  The  magnitude  of  the  ADC  error  can  be  checked  in  flight  by  observing 
the  winds  during  a 360  deg  level  turn.  If  the  actual  winds  do  not  change 
and  inertial  velocity  errors  do  not  change,  the  change  in  the  magnilude 
of  the  computed  winds  will  be  caused  by  an  ADC  error  (one  half  of  the 
observed  wind  change  is  the  ADC  error). 


Flat  Earth  Assumption 


55.  The  INS  gives  velocities  and  platform  altitudes  with  respect  to  a North 
pointing  local  vertical  (NPV)  axis  system.  This  axis  system  will  therefore  rotate 
and  experience  centrifugal  acceleration  as  the  result  of  the  earth's  rotation  and  the 
airplane's  translation  over  the  surface  of  the  earth.  The  INS  velocities  and 
altitudes  were  treated  as  inertial  quantities;  that  is,  these  quantities  were 
differentiated  to  obtain  acceleration  and  attitude  rates.  Strictly  this  is  only  valid 
if  these  quantities  were  obtained  with  respect  to  a nonrotating  nonaccelerating 
(i.e.,  fixed,  flat-earth)  axis  systems.  The  purpose  of  this  section  is  to  quantify  the 
effects  of  this  fixed  flat  earth  assumption  when  used  with  a North  pointing  local 
vertical  INS.  The  errors  introduced  by  making  a flat  nonrotating  earth  assumption 
when  using  a NPV  axis  system  can  be  visualized  by  the  following  discussion. 


56.  Round  Earth  - First,  it  is  assumed  that  the  earth  is  round  (spherical)  but  still 
nonrotating.  It  is  easy  to  visualize  the  following  two  effects  because  of  the 
translation  of  the  airplane  over  the  surface  of  the  earth. 


a.  Pitch  Rate  - Since  a slight  "pushover"  is  required  to  follow  the  curvature 
of  the  earth,  the  airplane  will  experience  a small  nose  down  pitch  rate. 
This  will  occur  at  a constant  indicated  aircraft  attitude  ( = 0)  with  respect 
to  the  NPV  axis  system. 
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b.  Vertical  Acceleration  - Because  of  this  slight  "pushover"  the  airplane  will 
experience  a small  acceleration  toward  the  center  of  the  earth  and  as  a 
result  the  load  factor  (N/’FP)  will  be  less  than  1 g.  This  effect  will  occur 
even  if  the  aircraft  mainieins  a constant  altitude  and,  therefore,  zero 
indicated  vertical  acceleration  (Vy)  with  respect  to  the  NPV  axis  system. 

57.  Round  Rotating  Earth  - Since  the  earth  is  rotating,  the  following  effects 
occur: 

a.  Additional  pitch,  roll,  and  yaw  rates  (U,  if>,  J1 ) caused  by  the  earth's  rotation 
rate  occur.  These  terms  exist  even  with  a constant  indicated  aircraft 
attitude  ( ,‘,l  , ) when  referenced  to  the  NPV  axis  system. 

b.  Centrifugal  acceleration  acting  orthogonal  to  the  spin  axis  of  the  earth 
occurs.  This  acceleration  lias  a North  and  Vertical  component;  however,  it 
would  not  be  evident  by  the  indicated  accelerations  from  the  NPV  axis 
system. 

c.  An  additional  Coriolis  term  is  caused  by  the  tangential  (or  Easterly) 
acceleration  of  the  airplane  because  of  the  change  in  radial  distance  from 
the  earth's  axis  of  rotation  (i.e.,  for  an  airplane  in  level  flight  the 
tangential  velocity  decreases  as  latitude  is  increased). 

58.  The  above  describes  the  errors  that  will  be  present  when  the  velocities  and 
attitudes  from  the  NPV  axis  system  are  treated  as  if  they  came  from  a 
nonaccelerating,  nonrotating  axis  system.  However,  the  magnitude  of  these  errors 
can  be  easily  calculated. 

Coriolis  and  Centrifugal  Acceleration  Effects 

59-  The  Coriolis  and  centrifugal  acceleration  effects  can  be  combined  to  produce 
equations  that  yield  the  magnitude  of  acceleration  errors  that  result  from  making 
the  flat,  nonrotating  earth  assumption.  The  equations  to  calculate  the  North,  East, 
and  Vertical  Acceleration  Errors  are: 

An  - Aj^'  = AA^  = 2 11  V^j  sin  (LAT) 

A£  - AE'  = AAe  = 2 Si  V^  cos  (LAT)  - 2 ,!  VNJ  sin  (LAT) 

Ay  - Ay  = ‘ Ay  = -2  Q V£].  cos  (LAT)  - . X (CX  R) 

j 

1 
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= North  Acceleration  Referenced  to  Inertial  Space 
= North  Acceleration  Referenced  to  the  NPV 


Ag  = East  Acceleration  Referenced  to  Inertial  Space 
Ag'  = East  Acceleration  Referenced  to  the  NPV 
Ay  = Vertical  Acceleration  Referenced  to  Inertial  Space 
Ay'  = Vertical  Acceleration  Referenced  to  the  NPV 
ft  = Earth's  Rotation  Rate  (7.29  x 10  ^ rad/sec) 

LAT  = Geodetic  Latitude 
X = Cross  Product 

60.  The  following  "trends"  can  be  derived  from  the  above  equations. 

a.  Most  of  the  error  terms  are  proportional  to  the  airplane's  inertial  velocity 
and,  in  general,  the  error  will  be  the  largest  for  an  airplane  at  high  speeds. 

b.  Latitude  (LAT)  and  the  airplane's  heading  (ip)  are  very  influential  factors. 

61.  The  following  examples  illustrate  the  errors  for  an  airplane  at  a true  airspeed 
of  661  KTAS  (1  224  km/hr)  or  Mach  1.0  at  sea  level. 

a.  North  Acceleration  Error  - The  largest  North  acceleration  errors  will 
occur  at  high  speeds  with  the  airplane  on  an  Easterly  heading  at  high 
latitudes.  For  example,  when  heading  East  at  Mach  1 at  40  deg  North 
latitude,  the  North  acceleration  error  would  be 

AAN  = .003  g 

b.  East  Acceleration  Error  - The  largest  East  acceleration  error  (AAg)  will 
occur  at  high  speeds  at 

Low  latitudes  with  the  airplane  in  a vertical  dive  ( y = 90°).  For 
example,  at  the  Equator  with  the  airplane  in  a Mach  1.0  vertical  dive, 
A A-  = .005  g. 


High  latitudes  with  the  airplane  heading  North.  For  example,  when 
heading  North  at  Mach  1 at  40  deg  North  latitude,  AAg  = .003  g. 

c.  Vertical  Acceleration  Error  - The  Vertical  Acceleration  Error  (AAy)  will 
be  largest  when  the  airplane  is  heading  East  at  low  latitudes.  For  example, 
at  the  Equator  when  heading  East  at  Mach  1.0,  AAy  = .0085  g. 
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62.  In  summary^  it  was  shown  that  errors  occur  when  the  fixed  flat  earth 
assumptions  are  applied  to  a NPV  referenced  axis  system.  For  a subsonic  airplane) 
the  maximum  error  introduced  will  be  .0085  g. 


63.  The  errors  discussed  in  this  section  were  considered  acceptable  for  this 
application;  however,  they  may  be  appreciable  for  other  applications)  e.g.,  high 
speed  at  critical  latitudes.  The  correction  terms  discussed  herein  can  be 
implemented  easily  in  the  software  program)  and  USNTPS  will  compensate  for 
these  errors  in  its  next  iteration  in  IFP  development. 

64.  The  following  block  diagram  (figure  8)  illustrates  how  the  acceleration 
corrections  (AAjj,  AA^,  A Ay)  can  be  converted  to  corresponding  body  referenced 

acceleration  corrections  (AA  , A A , A A ).  These  correction  terms  would  be  added 

x y z 

to  the  previously  computed  body  referenced  accelerations  (paragraphs  27  through 

31)  to  get  the  corrected  A , A , and  A . 

° x y z 


aan 

aae 


T A 

\ xi 

Fo 

V 


N AA, 


(Same  as  in 
paragraph  22) 


Figure  8 

Acceleration  Correction  Block  Diagram 
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Attitude  Rate  Effects 

65.  The  platform  attitudes  (^,  ij)  , 0 ) are  measured  with  respect  to  the  NPV  axis 
system.  As  previously  discussed,  this  axis  system  has  a rotation  rate  caused  by  the 
earth's  rotation  rate  and  the  aircraft's  translation  over  the  surface  of  the  earth. 
j The  differences  between  attitude  rates  ( 0,  ij  ) measured  with  respect  to  the  NPV 
axis  system  and  the  "true"  attitude  rates  are: 

AO-  -W  cos  (LAT)  sin  p - ± (VNJ2  + V^2)  1/2 
A<j>-  S2cos  (LAT)  cos  cos  0 + ft  sin  (LAT)  sin  0 
Ai p=  iJsin  (LAT) 


66.  The  magnitude  of  the  pitch  rate  error  is  illustrated  by  the  following 
example:  An  airplane  is  at  the  Equator  heading  East  at  600  KTAS  (949  km/hr)  with 
ii  = 0 . If  the  airplane  remains  at  the  above  conditions  (0  = 0),  the  attitude  rate 
(n),  as  indicated  by  the  NPV  data,  will  be  zero.  However,  the  inertially  referenced 
pitch  rate  will  not  be  zero.  The  inertially  referenced  pitch  rate  will  be  -.0042 
deg/s  (nose  down)  because  the  earth's  rotation  rate  plus  -.0028  deg/s  because  of  the 
aircraft's  translation  over  the  earth's  surface.  Thus,  a total  pitch  rate  error  of  .007 
deg/s  would  occur. 

67.  The  above  example  represents  a worst  case  condition  (i.e. , the  earth's 
rotational  term  and  the  aircraft's  translational  term  are  additive).  The  translation 
effect  will  be  greater  at  higher  true  airspeeds,  and  obviously  the  combination  of 
earth  rotation  and  aircraft  translation  will  be  lower  when  heading  West.  However, 
for  current  aircraft,  this  term  may  be  neglected. 

Differentiation  Errors 


68.  The  differentiation  of  the  INS  inertial  velocities  (Vj^j,  V^j,  V^j)  and  attitudes 

(i|M>  1 1|> ) was  required  to  get  the  IFP  data  (i.e.,  NZFP,  p,  q,  r,  etc).  In  general,  the 
errors  introduced  by  a smoothing/differentiation  routine  will  be  determined  by: 

a.  The  error  power  spectral  density  of  the  quantities  being  differentiated. 

b.  The  effective  bandwidth  produced  by  the  smoothing  routine  and  the  sample 
rate. 

It  is  usually  desirable  to  use  the  highest  possible  bandwidth;  however,  the  high 
frequency  noise  and  maximum  sample  rate  will  limit  the  bandwidth. 
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69.  The  "raw'1  inertial  velocities  and  attitudes  have  very  low  noise  levels,  and 
therefore  a relatively  high  bandwidth  smoothing  routine  was  used  (3  deg  polynomial 
through  5 dat. i points).  This  routine  gives  an  approximately  bandwidth  of  30 
rad/sec  with  the  sample  rate  of  25/sec  and  produces  low  noise  levels  (Appendix  C). 

70.  The  apparent  noise  level  exhibited  by  the  IFP  data  is  low.  The  noise,  however, 
appears  to  be  concentrated  at  high  frequencies  which  should  not  adversely  affect 
normal  airplane  stability  and  control  work.  Lower  noise  levels  may  be  obtained  by 
increasing  the  smoothing  {decreasing  the  bandwidth)  if  a lower  bandwidth  is 
permissible. 

ERROR  SENSITIVITY  ANALYSIS 


71.  The  previous  section  described  the  sources  of  errors  in  the  IFP  data.  The 
purpose  of  this  section  is  to  quantify  the  effects  of  various  INS  and  air  mass  errors 
on  the  IFP  data.  In  order  to  determine  how  an  error  affects  the  IFP  data,  one  can 
trace  the  error  through  the  various  transformations  that  were  used  to  calculate 
each  inertially  derived  parameter  (paragraphs  20  through  33).  Having  done  this, 
one  can  show  the  sensitivity  of  each  inertially  derived  flying  quality  and 
performance  parameter  to  each  error  source. 


72.  The  procedure  for  determining  the  sensitivity  of  an  inertially  derived 
parameter  to  an  error  source  is  in  concept  simple.  For  example,  to  determine  what 
angle  of  attack  error  (Au)  would  be  caused  per  unit  of  INS  velocity  error  (AVj),  one 
can  perturbate  the  inertial  velocity  (Vj)  by  the  assumed  error  (AVj)  and  then 
compute  the  resultant  change  in  the  computed  angle  of  attack  (a).  The  sensitivity 
term  is  then: 

<Ki  • A i 


The  angle  of  attack  error  is  therefore 


Au  = 


( 


tJ  K 


) A Vj 


73.  The  sensitivities  of  the  IFP  variables  to  the  primary  sources  of  errors  were 
computed  for  various  flight  conditions  (V^,  $ , 0 , i|/ ).  Table  II  summarizes  the 
results  of  this  sensitivity  analysis. 
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Table  H 

IFP  Error  Sensitivities 


IFP 

Variable 

Primary  Error 
Sources 

Sensitivity  (Note  2) 

Flight  Condition  to  Maximize 
Erro.  (i.e.,  Maximum  Sensitivity) 

< 

Erroneous  INS  Velocity 
Change  During  Maneuver 

. . 07.3  dec 

■V.  V.  kt 

I A 

1 

Total  INS  velocity  error  ( .’Vj)  is  j 
orthogonal  to  the  airplane's 
velocity  vector  and  the  airplane's 
y axis.  j 

Wind  Change  During 
Maneuver 

1 ‘ •_  57.3  deg  , 

•V  “ V kt 

W A 

Total  wind  vector  error  ( ' V^)  is 
orthogonal  to  the  airplane's 
velocity  vector  and  the  airplane's  | 
y axis. 

Sidpslip  Error  at 
Wind  Freeze 

= sin  ; cos  \ , 

in  deg/kt 

Airplane  remains  on  j^me  heading 
and  changes  bank  angle  by  90  ! 

after  wind  is  frozen  (Note  I). 

ADC  Error  at  Wind 
Freeze 

. 57.3  sin  . sin...  Airplane  changes  bank  n 

Y V ! angle  and  heading  hv  90 

in  deg/kt  A after  wind  is  frozen  (Note  1) 

» 

Erroneous  INS  Velocity 
Change  During  Maneuver 

1 57.3  deg 
ITT-  T~  Tt 

I A 

INS  velocity  error  vector 
( ' Vj)  occurs  along  airplane's 
v axis. 

Wind  Change  During 
Maneuver 

■ • 57.3  deg  , 

VW  VA  k‘ 

Total  wind  vector  error  ( ’V^l 
along  airplane's  y axis. 

i 

t 

Sideslip  Error  at 
Wind  Freeze 

. cos.. .cos'  ; 

“TT"  in  deg/deg 

Airplane  remains  at  same 
nominal  pitch  and  roll  angle 
after  wind  is  frozen  (Note  11. 

| 

ADC  Error  at  Wind 
Freeze 

~ = in  deg/kt  * 

' AC  ' A 

Airplane  makes  90°  heading 
change  after  freezing  winds 
(Note  1). 

AFPI 

1 

Angle  of  Attack  Error 

AFPI  . NZFP 
” 57.3 
in  g/deg 

Large  value  of  airplane 
acceleration  normal  to 
flight  path. 

| NZFP 

j 

. 

Angle  of  Attack  Error 

NZFP.  AFPI. 

, = 57.3  m g/deg 

_.  . . 

Large  valup  of  airplane 
longitudinal  acceleration. 

! p 

s 

i 

ADC  Error  at 
Wind  Freeze 

ps  • 1.69  AFPI 

V in  ft/sec/kt 

AC 

Large  value  of  airplane 
longitudinal  acceleration. 

Wind  Change  After 
Freezing  Wind 

Ps  . 1.69  NZFP  , 
(V’  in  ft/sec/kt 

Total  wind  vector  error 
orthogonal  to  flight  path 
and  airplane's  y axis. 

Erroneous  INS 
Velocity  Change 
During  Maneuver 

Ps  i 1.69  NZFP  t 
Vj  in  ft/sec/kt 

Short  term  INS  velocity 
vector  error  orthogonal  to 
flight  path  and  airplane’s 
y axis. 

1 

1 

Angle  of  Attack  Error 

P 1 AQ 

£_  - — N7FP  V 

- 57.3  VAC 

in  ft/sec/dea 

1 

Large  value  of  airplane 
acceleration  normal  to 
flight  path  and  high  true 
airspeed. 

* Worst  case  seasitivity,  i.e.,  the  error  vector  acting  orthogonal  to  or  along  a particular  axis. 

NOTES  I.  The  winds  are  frozen  with  ; ~ 0 and  , = 0. 

2*  The  units  of  all  velocities  (V^,  V^f  etc.)  should  be  in  knots. 
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74.  In  order  to  keep  table  II  to  a manageable  size,  it  does  not  present  the 
sensitivity  of  each  inertially  derived  parameter  to  every  error  source  but  rather 
only  the  significant  error  sources.  Many  inertially  derived  parameters  were 
omitted  from  table  13  since  their  error  sensitivities  are  obvious.  For  example,  the 
magnitude  of  A^,  A , and  A errors  will  have  a direct  correlation  (i.e.,  through  the 
transformations)  wi?h  the  magnitudes  of  A^,  Ag,  and  Ay  errors.  In  the  same  sense, 

the  magnitudes  of  V , V , and  V errors  have  a direct  correlation  with  the 
magnitudes  of  V^,  V^,  and  errors. 

75.  Table  II  shows  that  a and  3 errors  caused  by  erroneous  inertial  velocity 
changes  and  wind  errors  during  a maneuver  are  inversely  proportional  to  the  true 
airspeed.  For  example,  at  300  KTAS  (556  km/hr)  the  worst  case  u and  0 
sensitivities  to  short  term  INS  and  wind  errors  are  0.2  deg/kt  (.39  deg/m/s).  At  600 
KTAS  (1  111  km/hr)  the  worst  case  error  sensitivity  is  0.1  deg/kt  (0.19  deg/m/s). 


76.  Table  II  shows  that  P is  very  sensitive  to  short  term  inertial  velocity  errors, 
wind  vector  errors,  and  angle  of  attack  errors.  These  sensitivities  are  directly 
proportional  to  acceleration  normal  to  the  flight  path. 


! 


77.  An  analysis  of  the  error  sensitivities  shows  that  certain  maneuvers  can  be 
performed  which  will  maximize  an  inertially  derived  parameter's  error. 
Conversely,  maneuvers  can  be  selected  to  minimize  errors  in  the  parameters  of 
interest.  For  example,  with  an  ADC  error  the  angle  of  attack  errors  will 


be  minimized  by  conducting  maneuvers  on  approximately  the  same  heading  as 
existed  when  the  wind  was  frozen. 


FLYING  QUALITIES  DATA 


78.  Since  the  on-board  computer  sample  rate  of  all  inertial  parameters  was  25 
times  per  sec,  the  IFP  system  was  very  effective  for  stability  and  control 
testing.  At  this  sample  rate,  aircraft  frequency  response  up  to  30  rad/sec 
(approximately  five  cycles/sec)  could  be  adequately  documented.  This  frequency 
response  capability  should  be  adequate  for  most  aircraft  stability  and  control 
applications  and  is  well  within  the  20  cycle/sec  frequency  response  capability  of 
the  INS. 

79.  A random  sample  of  curves  showing  typical  stability  and  control  test 
maneuvers  is  presented  in  Appendix  C.  The  data  are  consistent  with  previously 
obtained  A-7  stability  and  control  data  and  the  cockpit  data  taken  during  the  test 
flights.  Of  note  is  that  each  of  these  curves  was  plotted  from  the  raw  unfiltered 
IFP  data.  As  previously  mentioned,  neither  cockpit  control  force  nor  control 
surface  position  data  were  available  during  this  feasibility  study  and  consequently 
are  not  depicted  on  any  of  the  curves  presented  herein. 
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80.  A data  sample  rate  of  fives  times  per  sec  appears  adequate  for  some  stability 
and  control  testing  such  as  dutch  roll  characteristics  (Appendix  C,  figures  1 thr<  'gh 
4).  However,  the  25  times  per  sec  sample  rate  is  required  for  documenting  the 
typically  higher  frequency  modes  of  motion  such  as  the  airplane's  short  period 
response  and  roll  performance.  Appendix  C,  figures  5 through  12,  give  comparisons 
of  the  same  short  period  maneuver  when  plotted  at  fives  times  and  25  times  per 
sec.  The  improved  quality  of  the  data  at  higher  sample  rate  is  apparent. 

81.  Acceleration  sensitivity,  or  the  variation  of  normal  acceleration  with  angle  of 
attack,  was  determined  by  plotting  acceleration  normal  to  the  flight  path  (NZFP) 
versus  ex.  Appendix  C,  figures  11,  12,  and  13,  show  NZFP  as  a function  of  a during 
a short  period  excitation  maneuver  and  during  a maneuvering  stability  wind-up 
turn. 

82.  Appendix  C,  figures  14  through  23,  present  time  histories  of  the  A-7C 
airplane's  roll  performance  at  200  KIAS  (103  m/s)  and  400  KIAS  (206  m/s).  The  :roll 
rates  and  the  magnitudes  and  direction  of  the  sideslip  during  these  roll  tests  were 
consistent  with  cockpit  data  and  available  A-7  stability  and  control  data. 

83.  A review  of  Appendix  C shows  that  all  parameters  are  essentially  noise  free. 
A comparison  of  Appendix  C,  figures  17,  18,  22,  and  23,  show  the  significant 
improvement  in  the  quality  of  the  angle  of  attack  measurement  when  the  inertially 
derived  a vice  the  airplane's  production  vane  a is  utilized.  The  noise  on  the  vane  a 
curves  appears  to  be  exacerbated  when  dynamic  pressure  and/or  angle  of  attack  is 
increased,  which  is  probably  a result  of  angle  of  attack  vane  flutter. 

PERFORMANCE  DATA 


84.  The  phase  of  flight  testing  to  which  IFP  provides  the  most  apparent  benefits  is 
aircraft  performance  testing.  A data  sample  rate  of  five  times  per  sec  was  used 
and  found  to  be  adequate  for  the  performance  testing  maneuvers  conducted  during 
this  program. 

TAKEOFF  AND  LANDING  PERFORMANCE 

85.  Takeoff  distances  can  be  determined  by  integrating  velocity  from  the  start  of 
the  takeoff  roll  until  weight-off-gear  or  until  vertical  velocity  starts  to  increase. 
Similarly  landing  distances  can  be  determined  by  integrating  velocity  from 
touchdown  (weight -on-gear  and/or  vertical  velocity  going  to  zero)  until  the 
airplane  stops.  If  the  inertial  system  has  residual  velocities  at  the  time  of  these 
tests,  the  data  can  easily  be  corrected  for  the  known  zero  velocity  which  exists  at 
the  beginning  of  the  takeoff  roll  and  at  the  end  of  the  landing  roll.  Using  the  IFP 
system  for  takeoff  and  landing  performance  tests  obviates  the  need  for  external 
camera  coverage  for  these  tests. 


26 


f, 


TM-TPS  76-1 


SPECIFIC  EXCESS  POWER 

86.  Specific  excess  power  (P  ) is  determined  with  conventional  instrumentation  by 

s 

taking  the  derivatives  of  pitot-static  derived  energy  height  (E^)  as  follows: 
h = altitude 

g = 32.2  ft/sec^  (9.81  m/sec^) 


v = true  airspeed  (ft/sec) 
y = flight  path  angle 


= h + 


dh  v dv 
dt  + g dt 


P dh  v dv 
S ' dt  + i dt- 


87.  Level  flight  acceleration  runs  are  typically  conducted  for  this  purpose.  With 
the  IFP  system,  Pg  can  be  determined  in  another  way. 

p - dh  y dv 

s df  g dt 

js  _ 1 dh  l dv 

V v dt  + g dt 

. . 1 dh 

and  since  — -tt-  = sin  Y 
v dt 


Q.V 

and  = AFP  (acceleration  along  the  flight  path) 


s . AFP 

— = sin  Y + 

v g 


and  sin  , + — = AFPI  as  defined  in  paragraph  30. 


Therefore,  P = V * AFPL 
’ s 


J 


r 

| ! 
f 
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88.  With  this  relationship,  P was  available  during  every  test  maneuver  by 
multiplying  indicated  acceleration  along  the  flight  path  (AFPI)  by  true  airspeed. 

89*  A level  flight  acceleration  run  was  conducted,  and  Appendix  C,  figure  24,  shows 
a comparison  of  Pg  derived  by  the  conventional  pitot-static  method  and  Pg  derived 
by  the  IFP  system.  The  data  show  good  agreement  over  the  entire  airspeed  range; 
however,  the  Pg  derived  by  the  IFP  system  appears  to  have  considerable  noise.  An 
examination  of  Appendix  C,  figure  25,  which  is  a curve  of  acceleration  normal  to 
the  flight  path  (NZFP)  versus  true  airspeed  (V^)  during  this  same  acceleration  run, 

shows  a high  degree  of  correlation  between  NZFP  and  the  P derived  by  the  IFP 

system.  That  is,  when  the  NZFP  shows  a peak  g greater  than  one,  the  Pg  curve 
shows  an  associated  "dip"  to  a local  minimum.  This  feature  of  the  data  is  exactly 
what  one  would  expect  (i.e.,  as  g increases,  the  airplane  drag  increases  and  thus  Pg 
decreases).  This  is  an  example  of  how  the  accuracy  and  dynamic  response 
capability  of  the  IFP  system  show  features  in  the  data  not  previously  discernible. 

90.  Appendix  C,  figure  26,  is  a plot  of  inertially  derived  angle  of  attack  versus 
vane  angle  of  attack  during  this  same  acceleration  run.  The  data  show  good 
correlation  at  the  low  angles  of  attack.  The  nominal  1.5  deg  difference  between 
the  angle  of  attack  values  from  each  source  is  probably  the  result  of  the  angle  of 
attack  vane  being  boresightea  to  a different  fuselage  reference  line.  No  further 
flights  were  conducted  to  confirm  this  hypothesis,  but  additional  testing  is 
recommended  to  develop  further  confidence  in  the  absolute  values  of  inertially 
derived  angles  of  attack. 

MANEUVERING  PERFORMANCE 

91.  A constant  airspeed  wind-up  turn  was  conducted,  and  the  variation  of  P as  a 

2 s 
function  of  (NZFP)  is  presented  in  Appendix  C,  figure  27.  As  expected,  the 

variation  of  P with  (NZFP)^  is  linear  at  the  low  g levels  and  then,  at  the  higher  g 
® 2 

levels,  decreases  more  rapidly  with  (NZFP)  . With  this  one  maneuver,  the  Pg 

characteristics  of  the  test  airplane  were  completely  documented  for  the  given 
altitude  and  airspeed.  From  these  data  the  following  were  determined: 

a.  The  level  flight  P (P  = 67.5  ft/sec  (20.6  m/s)  at  NZFP  = 1.0)  which  is 

s s 

usually  obtained  from  level  acceleration  runs. 

b.  The  sustained  g capability  of  the  airplane  (NZFP  = 2.59  g)  which  is  usually 
obtained  from  level  turn  performance  tests. 

c.  The  "high"  g Pg  values  (e.g.,  Pg  = 25,000  FPM  (127  m/s)  at  5 g's)  which  are 

not  available  from  conventional  flight  testing.  Conventional  methods  of 
determining  P at  different  g levels  require  knowledge  of  the  airplane's 
drag  polars. 
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DYNAMIC  PERFORMANCE  TESTING 

92.  The  purpose  of  dynamic  performance  testing  is  to  obtain  aircraft  performance 
information  from  dynamic  maneuvers  rather  than  from  time-consuming,  and 
therefore,  expensive,  steady-state  test  maneuvers.  It  has  been  estimated  that  a 70 
to  90  percent  savings  in  flight  time  could  be  realized  by  application  of  the  dynamic 
performance  method  in  an  extensive  performance  program  (reference  5). 

93.  The  basis  of  dynamic  performance  testing  is  the  accurate  determination  of 
flight  path  referenced  accelerations  (i.e.,  normal  to  and  along  the  airplane's  flight 
path).  With  conventional  dynamic  performance  methods,  special  acceleration 
packages  are  used  to  measure  the  accelerations  with  respect  to  the  airplane's  body 
axes.  These  accelerations  are  then  transformed  to  the  flight  path  axis  via  boom 
measured  angle  of  attack  (ctg)  and  sideslip  (3g).  These  boom  measured  quantities 

((Ttg,  3g)  are  subject  to  various  errors.  The  error  sources  and  the  causes  of  each  are 

listed  in  table  m. 


Table  HI 

Error  Sources  for  Boom  Measured  Angle  of  Attack  and  Sideslip 


Error  Source 

Causes 

Boom  Bending  Errors 

aE  = 

f (n_»  q > Q» 01  > M) 

BE  = 

f (AYI,  r,  Q,B,  M) 

Interference  Errors 

aE  = 

f (a,  3,  M) 

B E = 

f (a,  3 , M) 

Dynamic  Errors 

a E = 

f (q) 

ti 

w 

CQ 

f (r) 

94.  As  illustrated  in  table  m,  many  of  the  causes  of  errors  are  "dynamic"  in  nature 

(e.g.,  n , 0 , q,  r).  Approximate  corrections  are  applied  for  some  of  these  effects 
z 

(reference  5),  but  significant  residual  errors  still  occur.  Since  the  more  "dynamic" 
maneuvers  produce  larger  errors,  the  error  sources  become  a limiting  factor  for  the 
maneuvers  and  therefore  limit  the  potential  of  the  dynamic  performance  method. 
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95.  The  inertially  derived  a and  (J,  are  very  accurate  under  dynamic  conditions.  It 
was  shown  previously  that  the  primary  error  sources  for  the  IFP  data  jure  not  rate 
dependent.  In  fact,  some  of  the  primary  error  terms  are  minimized  if  the 
maneuver  is  completed  in  a short  period  of  time.  These  characteristics  of  the  INS 
derived  <t  and  (i  will  enhance  the  dynamic  performance  method  by  removing  some 
of  the  instrumentation  imposed  restrictions  on  the  flight  test  maneuvers.  This 
should  increase  the  efficiency  (amount  of  data  per  unit  flight  time)  of  the  dynamic 
performance  testing  method.  Additionally,  the  use  of  a proven,  highly  accurate  INS 
will  remove  the  need  for  the  extensive  calibration  required  with  a conventional 
dynamic  performance  package  which  sometimes  involves  calibrating  the 
accelerometers  before  and  after  each  data  flight. 

NEW  OR  IMPROVED  CAPABILITIES  FOR  FLIGHT  TESTING 


96.  Besides  providing  more  accurate  data  than  was  previously  obtainable,  IFP 
offers  additional  and  improved  flight  test  capabilities.  Many  of  these  new 
capabilities  obviate  the  need  for  external  data  gathering  sources  and  therefore 
provide  significant  cost  and  time  savings  to  the  flight  test  team. 

97.  As  described  previously,  takeoff  and  landing  distances  can  be  determined  by 
integrating  velocity  during  the  takeoff  roll  and  landing  roll.  This  eliminates  the 
need  for  external  camera  coverage  to  measure  takeoff  and  landing  distances. 

98.  Many  aircraft  carrier  suitability  tests  can  be  performed  using  solely  the  IFP 
instrumentation  system  as  the  data  source.  Altitude  lost  during  catapult  minimum 
end  airspeed  determinations  can  be  determined  by  integrating  vertical  velocity 
after  the  catapult  launch,  whereas  previously  the  "Mod  one  eyeball"  of  the  carrier 
suitability  engineer  when  combined  with  the  calibrated  "seat  of  the  pants"  feel  of 
the  test  pilot  was  utilized  to  determine  sink  off  the  bow.  (Historically,  camera 
coverage  produced  inaccurate  altitude  lost  results.) 

99.  In  carrier  suitability  testing  postcatapult  launch  longitudinal  acceleration, 
expressed  as  a/g  (longitudinal  acceleration  divided  by  gravity),  must  be  greater 
than  0.04  (reference  6).  Previously,  the  value  of  a/g  had  to  be  calculated  from 
thrust  required  and  thrust  available  data  (i.e.,  a/g  = (T-D)/W),  and  the  actual 
postcatapult  launch  a/g  was  never  directly  measured.  IFP  will  finally  give  the 
carrier  suitability  test  team  the  capability  to  measure  a/g  directly  even  during  the 
very  dynamic  catapult  launch  evolution. 

100.  The  altitude  lost  during  a wave-off  has  been  determined  traditionally  at 
NATC  by  using  SPN-42  radar  tracking.  With  IFP  the  same  degree  of  data  accuracy 
during  this  test  should  be  available  without  the  use  of  the  external  radar. 
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101.  Carrier  suitability  structural  landing  demonstrations  are  still  plagued  in  not 
being  able  to  measure  with  repeatable  accuracy  aircraft  landing  sink  speed  and 
aircraft  touchdown  attitude.  Camera  coverage,  doppler  radar,  and  TRODI 
(touchdown  rate  of  descent  indicator)  mirrors  have  all  been  tried,  but  none  have 
proven  to  have  the  desired  repeatable  accuracy.  IFP  offers  another  alternative 
which  should  prove  more  accurate  than  any  previous  method  and  it  requires  no 
external  tracking  device.  Since  IFP  lends  itself  well  to  telemetry,  real-time  data 
processing  could  obviously  be  used  for  these  tests. 

102.  Static  source  position  error  determinations  at  normal  accelerations  other 
than  one  and  during  transonic  flight  are  possible  from  using  the  IFP  system.  IFP 
will  not  replace  traditional  methods  for  airspeed  calibration  at  one  "g,"  but  IFP 
offers  a method  of  calibrating  static  source  position  error  during  acrelerated  flight 
which  was  never  previously  possible.  This  should  be  quite  valuable  in  the 
development  of  software  for  an  air-to-ground  weapon  system  that  uses  barometric 
inputs  in  its  weapons  delivery  computations. 

103.  As  discussed  in  the  Performance  Section,  IFP  is  an  outstanding  system  for 
dynamic  performance  testing.  Having  a continuous  readout  of  specific  excess 
power  (P  ) gives  the  flight  test  team  new  methods  to  obtain  performance  data  and 

enables  them  to  directly  obtain  performance  comparisons  that  previously  required 

the  knowledge  oi  the  aircraft's  drag  polars  (e.g.,  plotting  P at  different  g levels). 

s 

104.  A constant  Mach  wind-up  turn  or  a constant  g acceleration  run  can  be 
conducted  to  determine  Pg  as  a function  of  airspeed  and  normal  acceleration.  This 

will  obviate  the  need  for  having  the  test  pilot  conduct  the  traditional  stabilized 

point  level  turn  performance  tests.  (The  g at  which  P equals  zero  is  the  level 

s 

flight  sustained  turn  performance  at  each  specific  airspeed.) 

105.  Having  accurate  Pg  data  available  at  various  g levels  will  enable  the  flight 
test  team  to  better  compare  the  aircraft's  actual  performance  with  its  design 
criteria.  It  will  also  give  the  requirements  determination  people  the  opportunity  to 
measure  actual  maneuvering  performance  in  current  aircraft  against  potential 
threats  and  to  more  accurately  quantify  the  performance  needs  of  future  aircraft 
to  successfully  engage  potential  threats.  The  operational  requirements  people  will 
have  the  capability  of  more  completely  defining  a needed  performance  envelope, 
and  the  testers  will  be  better  able  to  document  each  aircraft's  total  performance 
envelope. 

106.  Since  IFP  continuously  calculates  flight  path  angle,  (y),  an  aircraft's  open 
loop  vertical  height  response  to  control  inputs  can  be  accurately  documented  (i.e., 
Bode  plots  can  be  made  of  longitudinal  control  input  frequency  to  output  gain). 
This  capability  should  permit  easier  derivation  of  the  software  to  be  used  in 
automatic  control  such  as  for  automatic  carrier  landings. 
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107.  Because  of  the  high  accuracy  and  low  noise  of  IFP  data,  it  should  be 
extremely  valuable  in  parameter  identification  work  (Note  1).  Flying  quality 
parameters  which  have  shown  themselves  as  being  noisy  in  their  application  to 
parameter  identification  are:  airspeed,  pitch,  roll,  and  yaw  accelerations;  normal 
acceleration;  and  longitudinal  acceleration.  Pitch,  roll,  and  yaw  accelerations  are 
particularly  difficult  to  obtain,  and  three  axes  angular  rate  measurement  devices 
arc*  currently  being  developed  for  this  purpose.  IFP,  however,  produces  the 
accurate,  low-noise  determination  of  all  the  necessary  aircraft  motion  parameters 
for  pa:ameter  identification  work.  It  is  recommended  that  a determination  of  IFP 
data  consistency  for  parameter  identification  work  be  conducted  such  as  by  using 
an  optimum  filtering/smoothing  routine. 


NOTF.  1.  Parameter  identification  or  parameter  estimation  are  the  methods  used 
in  system  identification  problems.  In  aircraft  flight  testing  this 
technique  has  traditionally  referred  to  the  determination  of  stability 
and  control  derivatives.  The  technique  involves  the  estimation  of  the 
characteristics  (e.g.,  stability  derivatives)  from  input/output 
measurements.  Usually  the  set  of  equations  that  describe  the  response 
being  observed  is  known,  and  parameter  identification  is  used  to 
estimate  the  value  of  each  term  in  the  set  of  equations. 

COST  COMPARISON 


108.  An  instrumentation  system  employing  the  IFP  concept  offers  a significant 
cost  savings,  especially  if  the  aircraft  being  tested  has  an  INS.  If  an  onboard  INS  is 
not  available,  the  relatively  high  cost  of  an  INS  may  dissuade  many  potential  users 
from  employing  IFP. 


109.  Since  rate  gyros,  accelerometers,  attitude  sensors,  a boom  with  angle  of 
attack  and  sideslip  vanes,  and  the  associated  signal  conditioning  is  not  required  if 
IFP  is  utilized,  the  attendant  cost  savings  are  obvious.  Since  these  items  do  not 
have  to  be  installed,  maintained,  or  calibrated,  there  is  also  a significant  savings  in 
the  labor  r ^uirements  to  install,  maintain,  and  calibrate  an  IFP  instrumentation 

system.  5 

j 

i 

110.  A nominal  ISO  man  days  of  labor  is  being  saved  in  installation  labor  at  TPS  by 

utilizing  IFP  in  an  A-7C.  This  number  was  obtained  by  comparing  the  IFP 
instrumentation  installation  labor  requirements  with  the  labor  requirements  for 
fully  instrumenting  a fixed  wing  airplane  for  FQ&P  testing  at  TPS  using 
conventional  instrumentation.  This  is  not  a very  good  comparison  since  the  IFP  i 

instrumentation  system  will  have  significant  additional  instrumentation  capabilities  • 

which  previously  have  not  been  utilized  at  TPS  (i.e.,  those  additional 
instrumentation  parameters  required  for  dynamic  performance  testing,  e.g.,  flight 

path  accelerometers). 
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111.  With  IFP  no  requirement  exists  to  know  the  relative  magnitudes  of  each  of 
the  parameters  to  be  recorded  prior  to  the  flight.  Therefore,  instrumentation 
preflight  time  is  reduced  sine**  the  instrumentation  parameter  sensitivities 
(expected  peak  to  peak  values  of  each  parameter)  do  not  have  to  be  sot  up 
according  to  the  type  of  flight  being  flown. 

112.  An  IFP  system  enables  the  test  aircraft  to  use  its  on-board  instrumentation 
system  to  determine  quantities  which  historically  have  been  determined  by  use  of 
external  sensors.  Examples  discussed  previously  include:  takeoff  and  landing 
distances,  wave-off  altitude  loss,  structural  landing  tests  requiring  touchdown 
attitudes  and  sink  speed,  altitude  lost  during  catapult  minimum  end  airspeed  tests 
.aboard  ship,  open  loop  aircraft  vertical  height  response  for  automatic  control 
testing,  and  transonic  static  system  position  error  determination.  Since  external 
cameras,  radar  tracking,  or  chase  aircraft  are  not  required  for  the  previously  listed 
tests,  the  cost  savings  thereby  accrued  are  a result  of  IFP  instrumentation. 


CONCLUSIONS 


GENERAL 

113.  An  inertial  system  can  be  successfully  used  to  derive  the  conventional 
parameters  which  describe  the  motion  of  an  aircraft.  This  method  of 

instrumentation  offers  significant  advantages  over  conventional  methods  of 
instrumentation. 

SPECIFIC 


114.  Traditional  flying  qualities  and  performance  data  were  easily  obtained  with 
higher  accuracy  and  less  noise  than  conventional  instrumentation  systems 
(paragraphs  78  through  89). 

115.  A boom  for  angle  of  attack  and  sideslip  vanes  was  not  required  to  obtain 
these  parameters.  The  accuracy  of  these  parameters  (especially  during  dynamic 
maneuvers)  from  the  IFP  system  was  greatly  improved  over  conventionally  obtained 
angle  of  attack  and  sideslip  (paragraphs  26  and  95). 

116.  If  an  inertial  system  is  already  available  in  the  test  aircraft,  a significant 
cost  savings  can  be  accrued  from  IFP  (paragraphs  108  through  112). 

1 17.  External  data  coverage  is  not  required  for  many  tests  that  previously  required 
camera  or  radar  tracking  (paragraphs  85,  100,  and  101). 

1 18.  New  test  capabilities  are  available  with  IFP  such  as  specific  excess  power 
determination  during  any  maneuver  and  static  source  position  error  determination 
at  normal  accelerations  other  than  one  (paragraphs  96  through  106). 
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119.  I'll'-  (l.il. i arc  of  tin*  quality  that  make  them  highly  useful  for  dynamic 
performance  testing  and  flying  quality  parameter  identification  work  (paragraphs 
95  and  107). 

120.  If  an  aircraft  has  an  INS  that  determines  North,  East,  and  Vertical  Velocities 

and  aircraft  attitude;  and  ADC  that  calculates  true  airspeed;  and  a recording 
device  for  these  parameters,  the  aircraft  processes  a full  IFP  package  which  is 
compatible  with  NATC  software  to  calculate:  <*  , > , p>  q,  r,  p,  q,  r,  A , A , A , 

n , p^,  etc.  (paragraphs  20  through  33).  x y z 

121.  Disadvantages  of  IFP  instrumentation  include: 

a.  Accurate  air  mass  referenced  airspeed  is  required  (paragraphs  21  and  47). 

b.  An  independent  method  of  determining  zero  sideslip  is  required 
(paragraphs  25  and  26). 

c.  If  the  wind  changes  during  a test  maneuver,  the  data  will  have  errors 
(paragraph  54). 

d.  Presently,  there  is  no  compensation  for  the  vertical  component  of  wind 
(paragraph  21). 

e.  High  cost  of  purchasing  an  INS  (paragraph  108). 

122.  Additional  advantages  of  IFP  instrumentation  include: 

a.  Shorter  instrumentation  installation  time  (paragraph  110). 

b.  Lower  instrumentation  installation  cost  (paragraph  1 09)- 

c.  Reduction  in  instrumentation  maintenance  (paragraph  109). 

d.  Calibrations  of  attitude,  acceleration,  and  rate  sensors  are  not  required 
(paragraphs  95  and  109). 

RECOMMENDATIONS 

123.  More  fully  utilize  the  capabilities  of  an  aircraft's  inertial  system  for  FQ&P 
testing. 

124.  Investigate  the  development  of  a portable  IFP  system. 

125.  Determine  IFP  data  consistency  for  parameter  identification  work 
(paragraphs  90  and  107). 
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TRANSFORMATION  OF  IMU  MOTION  TO  C.  G.  MOTION 


Vx  C.  G.  velocity  along  x axis 

Vy  C.  G.  velocity  along  y axis 

V C.  G.  velocity  along  z axis 

z 

V ' IMU  sensed  velocity  along  x axis 

Vy'  IMU  sensed  velocity  along  y axis 

V ' IMU  sensed  velocity  along  z axis 

z 

x distance  from  C.  G.  to  IMU  along  x axis  of  aircraft 

y distance  from  C.  G.  to  IMU  along  y axis  of  aircraft 

z distance  from  C.  G.  to  IMU  along  z axis  of  aircraft 
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P«  q»  r 


AV  , AV  , AV 
x y 2 


AA  , A A ,AA 
x y z 


A V 


AV 

A V_ 


A A 


A A 


A A 


true  angular  rates  of  aircraft  about  C.  G. 

im  ~ rnents  in  sehsed  body  axes  velocities  at  the  IMU 
because  of  aircraft  rotation 

increments  in  sensed  body  axes  accelerations  at  the 
IMU  because  of  aircraft  rotation 

y r - z q 

z p - x r 

x q - y p 

2 2 

yr-zq  + x(q  +r)-p(yq  + zr) 

2 2 

z p - x r + y (p  + r ) - q (z  r + x p) 

2 2 

x q - y p t-  z (p  + q ) - r (x  p + y q) 


Assuming  no  bending  and  a rigidly  mounted  IMU,  the  attitude  rates  and 

attitude*  accelerations  of  the  IMU  will  be  identical  to  those  of  the  C.  G.  regardless 

of  IMU  location  on  the  aircraft.  If  the  INS  gives  a velocity  output,  accelerations 

are  determined  by  differentiating  velocity  along  each  axis.  Therefore,  only  V , V , 

and  V need  to  be  corrected  if  the  IMU  and  C.  G.  are  not  collocated, 
z 

In  order  to  correct  the  IFP  derived  velocities  about  the  IMU  to  velocities  about 
the  C.  G.,  the  increments  in  body  axes  velocities,  calculated  above,  must  be  added 
to  the  IMU  sensed  velocities. 


V = V ’ + AV 

XX  X 


V 

y 


V 1 + AV 

y 


y 


V = V ’ + AV 
z z z 


Accelerations  about  the  IMU  could  be  corrected  in  a similar  way,  but  as 
explained  previously,  this  isn't  required.  These  same  equations  obviously  can  be 
used  to  determine  the  velocities  and  accelerations  about  any  known  point  on  the 
aircraft. 
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FLYING  QUALITY  AND  PERFORMANCE  CURVES 


ROLL  RATE 
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MODEL  A-7C  AIRPLANE 
BuNo  156776 


-60.00  •— 
0.00 


CONFIGURATION 
DATA  SAMPLE  RATE 
ALTITUDE 
AIRSPEED 
YAW  STAB 


CRUISE 

5/SEC 

20,000  FT  HP 
200  KIAS 
OFF 


6.00 

TIME  - SEC 


METRIC  CONVERSION 

1000  ft  « 304.8  m 
100  kt  - 51.44  m/s 


Figure  1 

Dutch  Roll  Excitation 
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Figure  2 

Dutch  Roll  Excitation 
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ANGLE  OF  ATTACK  (INS) 
XALPHA 
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MODEL  A-7C  AIRPLANE 
BuNo  156776 


h 

f t 


CONFIGURATION:  CRUISE 

ALTITUDE:  20,000  FT  HP 
AIRSPEED:  300  KIA3 
CONTAUG:  OFF 
DATA  SAMPLE  RATE:  25/SEC 


\J 


6.00 

TIME  • SEC 


METRIC  CONVERSION 

1000  FT  = 304.8  M 
100  KT  = 51.44  M/S 


Figure  5 

Short  Period  Excitation 


43 


APPENDIX  C 


ANGLE  FOR  ATTACK  (INS) 
XALPHA 
DEG 
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MODEL  A-7C  AIRPLANE 
BuNo  156776 


CONFIGURATION:  CRUISE 

ALTITUDE:  20,000  FT  HP 
AIRSPEED:  300  KIAS 
CONTAUG:  OFF 
DATA  SAMPLE  RATE:  5/SEC 


V 


6.00 

TIME  • SEC 


METRIC  CONVERSION 

1000  ft  = 304.8  m 
100  kt  = 51.44  m/s 


Figure  6 

Short  Period  Excitation 
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MODEL  A-7C  AIRPLANE 
BuNo  156776 


CONFIGURATION:  CRUISE 

ALTITUDE:  20,000  FT  HP 
AIRSPEED:  300  KIAS 
CONTAUG:  OFF 
DATA  SAMPLE  RATE:  25/SEC 
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0.00 


6.00 

TIME  - SEC 


METRIC  CONVERSION 

1000  ft  = 304.8  m 
100  kt  = 51.44  m/s 


Figure  7 

Short  Period  Excitation 
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PITCH  RATE  - 
DEG/SEC 
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Figure  8 

Short  Period  Excitation 


46 


APPENDIX  C 


ACCELERATION  NORMAL  TO  FLIGHT  PATH 
NZFP 


TM-TPS  76-1 


MODEL  A-7C  AIRPLANE 
BuNo  156776 


1.751 


1.50 


1.25 


O 1.00 


0.75 


0.501 


0.25 


0.00 


CONFIGURATION: 
ALTITUDE: 
AIRSPEED: 
CONT  AUG  . 
DATA  SAMPLE  RATE: 


CRUISE 
20,000  FT  HP 
300  KIAS 
OFF 
25/SEC 


2.00 


4.00 


6.00 

TIME  • SEC 


8.00 


10.00 


12.00 


METRIC  CONVERSION 

1000  ft  = 304.8  m 
100  kt  = 51.44  m/s 


Figure  9 

Short  Period  Excitation 
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Figure  10 

Short  Period  Excitation 
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Short  Period  Excitation 
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Roll  Performance 
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Roll  Performance 
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Roll  Performance 
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Ren!  Performance 
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Figure  23 
Roll  Performance 
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Specific  Excess  Power 
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Angle  of  Attack  Correlation 
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Maneuvering  Performance 
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TM-TPS  76-1 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 

ADC  - Air  Data  Computer 

AFP  - Acceleration  Along  the  Flight  Path 

AFPI  - Acceleration  Along  the  Flight  Path  - Indicated 

AYI  - Indicated  Lateral  Acceleration 

A - Acceleration 

A^  - Acceleration  Along  the  North  Axis 

A_.  - Acceleration  Along  the  East  Axis 

Ay  - Acceleration  Along  the  Vertical  Axis  (positive  down) 

A - Acceleration  Along  the  Airplane's  x Axis 

Ay  - Acceleration  Along  the  Airplane's  y Axis 

A^  - Acceleration  Along  the  Airplane's  z Axis 

CEP  - Circular  Error  Probable 

D - Airplane  Drag 

FQ&P  - Flying  Qualities  and  Performance 

2 2 

g - Gravitational  Constant  (32.2  ft/sec  or  9.81  m/sec  ) 

hr  - Hour 

H - Pressure  Altitude 

P 

IFP  - Inertially  Derived  Flying  Qualities  and  Performance  Parameters 
IMU  - Inertial  Measurement  Unit 

INS  - Inertial  Navigation  System 

kg  - Kilograms 

km  - Kilometer 
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kt 

KTAS 

M 


>n 


N 

NATC  - 

nmi 

NPV 

NZFP  - 

P 

P 

s 

P 

Q 

q 

R 


s 


s 


T 


USNTPS  - 
V 


V 


NI 


VEI 


VI 


VT 


Knot 

True  Airspeed  - Knots 
Mark  Number 
Meters 

Acceleration  Measured  in  "g"  Units 
Naval  Air  Test  Center 
Nautical  Miles 

North  Pointing  - Local  Vertical  Axis  System 
Acceleration  Normal  to  the  Flight  Path  (in  "g"  Units) 
Position 

Specific  Excess  Power 

Roll  Rate  Measured  About  the  Airplane's  x Axis 
Dynamic  Pressure 

Pitch  Rate  Measured  About  the  Airplane's  y Axis 
Earth's  Radius 

Yaw  Rate  Measured  About  the  Airplane's  z Axis 

Laplace  Operator 

Seconds 

Thrust 

U.  S.  Naval  Test  Pilot  School 
Velocity 

Inertial  Velocity  - North 
Inertial  Velocity  - East 
Inertial  Velocity  - Vertical 
Airplane's  Inertial  Velocity  Vector 
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TM-TPS  76-1 


V 


N 


V. 


V 


A 


AC 


V 


WN 


V 


WE 


WV 


V 


w 


V 

x 

V 

y 

v 

z 

W 


x 

y 

z 


- Air  Mass  Referenced  Velocity  - North 

- Air  Mass  Referenced  Velocity  - East 

- Air  Mass  Referenced  Velocity  - Vertical 

- Magnitude  of  Air  Mass  Referenced  Velocity  (True  Airspeed) 

- Airplane's  Air  Mass  Referenced  Velocity  Vector 

- Magnitude  of  Air  Mass  Referenced  Velocity  - Computed  from  "frozen" 
winds  and  inertial  velocities 

- Wind  Velcoity  - North 

- Wind  Velocity  - East 

- Wind  Velocity  - Vertical 

- Wind  Velocity  Vector 

- Component  of  Air  Mass  Referenced  Velocity  Along  Airplane's  x Axis 

- Component  of  Air  Mass  Referenced  Velocity  Along  Airplane's  y Axis 

- Component  of  Air  Mass  Referenced  Velocity  Along  Airplane's  z Axis 

- Airplane  Gross  Weight 

- Airplane's  x Axis 

- Airplane's  y axis 

- Airplane's  z Axis 


‘I 

‘B 

'e 


- Angle  of  Attack 

- Inertially  Derived  Angle  of  Attack 

- Boom  Measured  Angle  of  Attack 

- Angle  of  Attack  Error 

- Sideslip  Angle 

- Inertially  Derived  Sideslip 
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„ - Boom  Measured  Sideslip 

O 

, j,  - Sideslip  Error 

. - Flight  Path  Angle  (Air  Mass  Referenced) 

, - Flight  Path  Angle  (Earth  Referenced) 

, - Air  Density 

Special  Symbols 

(’ ) - Indicates  First  Derivative  of  Term  in  Parenthesis 

("  ) - Indicates  Second  Derivative  of  Term  in  Parentheses 

A - Change  or  Error  in 

( ) - Indicates  Vector 

( ) - Indicates  a Measured  Quantity 

( ) - Indicates  Value  of  Quantity  at  "Wind  Freeze"  or  Time  Zero 

X - Cross  Product 
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